Correlations in intermediate-energy two-proton removal reactions 
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We report final-state-exclusive measurements of the light charged fragments in coincidence with 
26 Ne residual nuclei following the direct two-proton removal from a neutron-rich 28 Mg secondary 
beam. A Dalitz-plot analysis and comparisons with simulations show that a majority of the triple- 
coincidence events with two protons display phase-space correlations consistent with the (two-body) 
kinematics of a spatially-correlated pair-removal mechanism. The fraction of such correlated events, 
56(12) %, is consistent with the fraction of the calculated cross section, 64 %, arising from spin 
S — two-proton configurations in the entrance-channel (shell-model) 28 Mg ground state wave 
function. This result promises access to an additional and more specific probe of the spin and 
spatial correlations of valence nucleon pairs in exotic nuclei produced as fast secondary beams. 

PACS numbers: 24.10.-i 24.50.+g 25.60.Gc 29.38.-c 



Access to nuclear reactions that can probe the states of 
pairs of nucleons in an atomic nucleus is a long-standing 
ambition. Specifically, an ability to probe the spin- 
structure of nucleon pairs and, for example, to identify 
and quantify two- nucleon correlations in the [S, T] = [0, 1] 
spin-isospin channel, is required. Intermediate-energy re- 
actions that remove two nucleons (2N) suddenly from a 
fast projectile, in collisions with a light target nucleus, 
have been shown [l|-[4| to provide a sensitivity to the 2N 
configurations near the projectile surface. Experimen- 
tally, the importance of nucleon removal reactions derives 
from their high detection efficiency (forward focusing of 
the reaction residues) and use of thick reaction targets, 
increasing the effective luminosity of the relatively low 
intensity exotic beams. Prior to the present work these 
2N removal cross section measurements have been inclu- 
sive with respect to the final states of both the target and 
the removed nucleons but often exclusive with respect to 
the bound final-states of the forward-traveling projectile- 
like residues, allowing spectroscopic studies based on the 
different residue final-stateyields and the shapes of their 
momentum distributions [5l4l0j. 

The present work sacrifices, temporarily, this better- 
understood residue final-state sensitivity to investigate 
new information that might be forthcoming from more 
exclusive measurements of the final states of the removed 
nucleons. The 9 Be( 28 Mg, 26 Ne) reaction was used at an 
intermediate energy of 93 MeV/u. The chosen reaction 
was studied previously in a 26 Ne-y coincidence measure- 
ment [Iol | and was used to confirm the predicted relative 



populations of the four bound 26 Ne final states [lj]. The 
present data set has also been used ll| to confirm that 
the measured contributions to the inclusive 2N-removal 
cross section from each of the possible clastic and inelas- 
tic removal mechanisms [l[ were consistent with calcu- 
lations that use eikonal reaction dynamics and sd-shell- 
model structure inputs for the 28 Mg to 26 Ne( J 7r ) 2p over- 
lap functions. There it was shown that only 8(2) % of 
the inclusive 9 Be( 28 Mg, 26 Ne) reaction cross section re- 
sults from elastic breakup, the remainder being associ- 
ated with reactions in which at least one of the removed 
protons interacts inclastically with the 9 Be target nu- 
cleus [m. 

The experiment was carried out at the Coupled Cy- 
clotron Facility at the NSCL(MSU). The 28 Mg secondary 
beam was produced by projectile fragmentation of a 140 
MeV/u 40 Ar primar y be am and selected using the A1900 
fragment separator [l2j. The 9 Be reaction target had a 
thickness of 100 mg/cm 2 and was placed at the target 
position of the high-resolution S800 magnetic spectro- 
graph The beam energy at mid-target position was 
93 MeV/u. The 26 Ne reaction residues were identified 
by measuring the energy loss in the ionization chamber 
in the S800 focal plane and the time of flight measured 
between scintillators before and after the reaction target. 
The 26 Ne energy and momenta were reconstructed from 
the measured positions and angles in the focal plane us- 
ing the position-sensitive cathode readout drift chambers 
(CRDC) and ion optical ray-tracing. Protons and other 
light charged particles were detected and identified in the 
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high-resolution array HiRA [14| . Other experimental dc- ^ 



tails and the angular coverage of HiRA were discussed in 



detail in Ref. [Ill ], In the present analysis we considered 
those triple-coincidence events where both of the light 
particles were identified as protons. These 4810 events 
represent 28(3) % of the inclusive two-proton removal 
cross section |ll| . 

The continuum energy of these 26 Ne+p+p triple- 
coincidence events in their total momentum P c \2 = 
frame can be reconstructed from the three measured lab- 
oratory frame four-momenta P\. Thus, 



(1) 



is the energy (above threshold) of the dissociated 28 Mg 
fragments, whose distribution is shown in Fig.[l]a). This 
distribution served as an input for Monte-Carlo kinemat- 
ics simulations, used to populate the 26 Ne+p+p phase 
space and explore correlations of the detected protons. 

In the simulations the three particles were assumed 
to be structure- less. There was no consideration of spin 
degrees of freedom of the two protons. Two kinematics 
scenarios were used, (a) three-body and (b) two-body, to 
generate the four-momenta of the three particles in the 
Pen — frame. The summed energies of the three par- 
ticles were sampled from the experimental da/dE b s dis- 
tribution shown in Fig. Ufa). In the three-body case the 
available energy -E bs was distributed among the three 
particles (democratic breakup). The residue and pro- 
tons are then correlated only by energy and momentum 
conservation leading to the broad two-proton relative en- 
ergy distribution shown in Fig. [TJb) . The two-body case 
simulated events where the target was assumed to have 
encountered spatially-localized proton pairs (carrying en- 
ergy e* in their P\i = frame) and delivered an impulse 
to these pairs; as would be expected for some fraction of 
events in the 2N-removal reaction mechanism that favors 
the surface localization and spatial proximity of pairs of 
nucleons (see e.g. [H). The energy £? ODS was now shared 
between an assumed e* and the motion of the center- 
of-mass of the protons, producing 26 Ne+(2p)* configura- 
tions that decay to 26 Ne+p+p. 

To study the two-proton correlations we first recon- 
struct the measured two-proton relative energy, E re \, dis- 
tribution shown in Fig. [TJb). These data points were 
fitted to extract the relative numbers of three-body (red) 
and two-body (green) kinematics events, the events for 
E Te ] > 40 MeV being assumed to result from the three- 
body mechanism only. The fitted two-body events then 
determine the e* spectrum that is sampled to generate 
the two-body phase-space. We attribute the fraction 
0.56(12) to these latter, correlated 2p-removal events. 
The two kinematics scenarios also show a different sig- 
nature with respect to the angle between the protons. 
However, due to both the coarse and somewhat limited 
azimuthal angles coverage in the experiment, no addi- 
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FIG. 1: (color online) (a) Energy of the dissociated 26 Ne+p+p 
fragments, in their P c \2 = frame. The cross section has been 
corrected for the geometric efficiency of the HiRA array, as 
described in [llT |. (b) Relative energy distribution of the two 
detected protons. The fractions of simulated three-body (red) 
and two-body (green) kinematics events have been fitted to 
the data set (see text for details). The inset compares the 
measured relative energy distribution with that derived from 
mixing protons from two independent events (red). The dis- 
tribution obtained from mixed events has been scaled to the 
experimental two-proton relative energy distribution above 
40 MeV 



tional useful evidence could be extracted from the mea- 
sured angular distributions. We note that an E ro \ distri- 
bution generated from mixed events - in which the two 
protons were completely uncorrelated - did not show the 
increase in the cross section for small E re \ , as is shown in 
the inset to Fig. Hfb). 

Sequential two-proton removal events, via a single- 
proton removal to proton-unbound intermediate-states 
in 27 Na, are not considered to contribute. Such indi- 
rect paths would involve the proton decay of intermedi- 
ate 27 Na excited states with energies in excess of 13.3 
MeV, whereas the neutron threshold in 27 Na lies below 7 
MeV in excitation energy [l(J. This expectation is con- 
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firmed experimentally by the analysis of the Dalitz plot 
in Fig. [2j as will now be discussed. 

In high-energy physics, three-particle final-states are 
commonly studied by the analysis of Dalitz plots flol ]. 
This shows the correlation of the squares of the invariant 
masses My = (Pi+Pj) 2 . In the absence of sources of two- 
particle correlations, the observed distributions will be 
uniform within the boundaries defined by energy and mo- 
mentum conservation. Since the quantity on the Dalitz 
plot is proportional to the square of the matrix element 
for the system, intermediate states as well as symmetries 
generate a non-uniform distribution. In the present case, 
the boundaries will be different for each value of the con- 
tinuum energy, E \> s , so we have adopted the approach of 
Ref. M and used the normalized invariant masses, Wy, 
that range from to 1, i.e. 



W 2 = 



M 2 



(rrii + rrij) z 



. (2) 

(i?obs + m% + mj) 2 - (m.( + m 3 ) 2 

The Dalitz plot of the current 26 Nc+p+p event data 
is shown in Fig. [H a). The experimental results are com- 
pared to the three-body and two-body 2p-removal kine- 
matics simulations using a two-body fraction of 0.56 in 
panel (d). The three-body case results in a flat distribu- 
tion while the two-body kinematics simulation produces a 
non-uniform filling of the phase space - with an increased 
intensity at lower proton-proton invariant masses Wp-p. 

If present, sequential two-proton removal via an ex- 
cited intermediate state in 27 Na would lead to vertical 
bands with constant W 2 p , as is shown in panel (e). This 
was not observed in the experimental data, providing 
strong support for our assumption that such indirect 
paths do not contribute. This provides a direct exper- 
imental confirmation of earlier arguments, based on en- 
ergetics, the cross section, and the inclusive momentum 
distribution, for the dominance of the direct removal of 
the two strongly-bound valence particles 0, ■ 

Also shown, in panels (b) and (c) of Fig. [2j are the 
data projections onto the W 2 p and W 2 p axes. These 
data projections are also reasonably well described by 
the combination of the two breakup modes determined 
from the fit to the relative energy spectrum, Fig. [TJb). 
We note that the effect of the 9 Be reaction target enters 
the simulations only very indirectly, through the mea- 
sured projectile excitation energy distribution, dcr/dE b s 
of Fig. [TJa). Thus, the comparisons of the data with the 
two simulated kinematics scenarios does not provide in- 
formation on the final states of the target, to which our 
results remain inclusive. 

It is instructive to compare the deduced correlated 2p- 
fraction to the calculated two-nucleon spin contributions 
to the reaction cross sections. The partial cross sections 
to each 26 Ne( J w ) final state receive incoherent (additive) 
contributions from the total orbital L and spin S angu- 
lar momentum components in the ( 26 Ne( J 7r )| 28 Mg) two- 
proton overlaps 01 . Table U shows the S=0 fractions (%) 
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FIG. 2: (color online) Dalitz plots of core ( 26 Ne residue)- 
proton W cp and proton-proton W pp invariant masses. Panel 
(a) shows the experimental data with projections on Wp P (b) 
and Wcp (c). The results of the simulation using a two- 
body fraction of 0.56 are shown in panel (d). Panel (e) 
shows the corresponding Dalitz plot for simulated sequential 
two-proton removal events, via a single-proton removal to a 
proton-unbound intermediate-state in 27 Na. 



of both the overlaps and the removal cross sections, as=o, 
when using the sd-shell-model two-nucleon amplitudes 
tabulated in Ref. [17|. These [S, T] = [0, 1] wave function 
components will have significant but not exclusively sq 
2p configurations. 

The relative transparency of the direct reaction mech- 
anism to the S of the two nucleons delivered by the pro- 
jectile is evident from Table HI The S = fractions of 
the overlaps are reflected rather directly, with only small 
final-state-dependent enhancements, in their percentage 
contributions to the cross sections, <rs=o ■ Here the S=0 
terms are seen to be responsible for 64 % of the com- 
puted inclusive cross section, consistent with the two- 
body event fraction deduced from the simulations. This 
suggests that the [S, T] = [0, 1] correlations present in that 
part of the 28 Mg ground-state wave function sampled by 



TABLE I: Calculated partial and inclusive cross sections for 
the 28 Mg(— 2p) reaction at 93 MeV per nucleon and the per- 
centage contributions of 5=0 terms to the two-proton over- 
laps and the cross sections. The 2p stripping contributions are 
shown in Table III of Ref. [13], there for 82.3MeV/nucleon. 
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the target may be maintained in the sudden two-proton 
removal process and be reflected, here via the two-body 
kinematics events simulation, as two-proton correlations 
in the final state. 

This being the case, an expectation is that the width of 
the parallel momentum distribution of residues, gated on 
the two-body- like events with smaller E re \, should be nar- 
rower; this due to the size of the 26 Ne(0 + ) ground state 
cross section with its narrow momentum distribution Q . 
However, as all 26 Ne final states contribute significantly 
to the 64 % 5 = contribution to the inclusive cross sec- 
tion, this differential width is only slightly larger than our 
current experimental resolution. Fig.[3]shows the parallel 
momentum distribution of the 26 Ne residues (in coinci- 
dence with two protons) gated on two-body and three- 
body events. We only consider the central part of mo- 
mentum distributions since the high and low-momentum 
tails were influenced by acceptance cut-offs and proton 
detection thresholds. The parallel momentum distribu- 
tion gated on three-body mechanism events (red line in 
Fig. |3]) is slightly wider than the corresponding distribu- 
tion gated on two-body events. This differential width 
effect is consistent with what is expected based on the 
theoretical 5 = and 5 = 1 momentum distributions 
folded with the width of the momentum distribution of 
the unreacted 28 Mg beam. Since the two- and three-body 
simulations describe the available data, and the extracted 
momentum distributions are consistent, within statistics, 
with the two-body events being closely allied with the 
5 = overlap function components, we may expect that 
the less-spatially-localized 5 = 1 two-proton components 
will track the results of the three-body simulation. Addi- 
tional data, e.g. to states with different 5 = 0, 1 fractions 
will be needed to confirm this expectation. 

In summary, we have observed significant kinematical 
correlations of two final-state protons measured in coin- 
cidence with 26 Ne residues following the 9 Be( 28 Mg, 26 Ne) 
two-proton removal reaction. We attribute the deduced 
56(12) % fraction of two-body events to the observed 
cross section as due to spin 5 = two-proton configura- 
tions in the 28 Mg wave function in the entrance-channel. 
This result suggests the potential for such measurements 
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FIG. 3: (color online) Parallel momentum distribution of the 
26 Ne residue gated on two-body (black) and three-body (red) 
events. The two components were extracted by applying a cut 
on E IB i at the intersection of the two distributions in Fig. [Tib). 
Contributions from three-body events at low E le i were sub- 
tracted from the two-body distribution and vice-versa for the 
two-body events with high relative energy. The lines are 
Gaussian fits. The two distributions have been normalized to 
show the difference in width. The blue distribution illustrates 
the width of the momentum distribution of the incoming 2S Mg 
beam. 



to provide an additional, more specific probe of the spin 
correlations of valence nucleon pairs in exotic nuclei. It 
also suggests the need for additional exclusive measure- 
ments to confirm and quantify this proposed spin sen- 
sitivity. Specifically, two-nucleon removal reactions to 
nuclei with only a single (0 + ) bound state, where 5 = 
configurations will dominate, would be ideally suited for 
this kind of study. Further indirect indications of 5 = 
driven final-state 2p correlations could come from obser- 
vations of enhanced 3 He and a-particle yields produced 
following In- or 2n-pickup from the target by 5 = 2p 
configurations in their surface-grazing removal collisions. 
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